Efimov states near a Feshbach resonance 
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We describe three-body collisions close to a Feshbach resonance by taking into account two- 
body scattering processes involving both the open and the closed channel. We extract the atom- 
dimer scattering length and the three-body recombination rate, predicting the existence at negative 
scattering length of a sharp minimum in the recombination losses due to the presence of a shallow 
bound level. We obtain very good agreement with the experimental results in atomic ^''"^Cs of 
Kraemer et al. [Nature 440, 315 (2006)], and predict the position of Efimov resonances in a gas of 
^^K atoms. 
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I. INTRODUCTION 

In the early seventies Efimov predicted that bosons 
with strong two-body interactions may bind together to 
form trimers with interesting universal properties, even 
under conditions where two-body bound states are not 
allowed In particular, an infinite series of weakly 
bound three-body levels appears close to the continuum 
threshold when the magnitude of the atom-atom scatter- 
ing length \a\ becomes large compared to the character- 
istic range of the potential i?,, which is of order lOOao for 
alkali atoms (ap being the Bohr radius). In the region 
close to a resonance of a, these states are universal in the 
sense that they map onto each other via a discrete scaling 
symmetry, obtained by multiplying lengths and energies 
by a factor A and 1/A^ respectively, where A ~ 22.7 
The observation of Efimov physics has been a major chal- 
lenge in modern physics, due to the smallness of the typi- 
cal binding energies. Experimental evidence of an Efimov 
state has been recently reported by Kraemer et al. in a 
thermal gas of ^'^^Cs atoms 3]. On the theoretical side, a 
microscopic description of Efimov physics is known to be 
particularly challenging. Indeed, the two-body scattering 
length a alone is not sufficient to determine the solution 
of the three-body bosonic problem. Remarkable results 
have been obtained within the so-called "universal the- 
ory" (UT) 0, which assumes \a\ ^ R and a two-body T- 
matrix proportional to {l + ika)~-^ , where i? — h^k^/m is 
the energy in the center-of-mass frame of the two atoms. 
Nonetheless UT has a priori no predictive power, since 
it contains two undetermined parameters, the first fix- 
ing the structure of the energy levels, and the second 
describing recombination into deeply-bound levels. 

The analysis of current experiments in ultracold gases 
by means of UT might not properly work for different rea- 
sons. In first instance, experiments are usually performed 
outside the regime \a\ ^ i?, in order to avoid large three- 
body losses which in bosonic systems rise rapidly as 
close to a resonance. In addition, most of the interest- 
ing Feshbach resonances are characterized by background 
scattering lengths |abg| which are large compared to the 
range R, signaling the presence of a shallow (long-lived 
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FIG. 1: (Color online) Two and three-body energy levels as a 
function of magnetic field for ^''"^Cs. In order to visualize the 
Efimov states, both axes display highly nonlinear scales, ob- 
tained by plotting {B-BqY^^ and E^^^° . The cross shows the 
position of the recombination resonance observed in Ref. 



or bound) state in the open channel of the resonance 
which is not taken into account by the scattering length 
approximation. It is therefore interesting to develop a 
model that does not rely on the assumptions of UT and 
compare it with the existing theories, such as those based 
on the hyperspherical method 0, Q or the AGS ap- 
proach [7|,|8|. 

We choose here to work in momentum space, where 
the energy dependence of two-body interactions can be 
included in a natural way. We incorporate low-energy 
collisions between three identical bosons by means of 
a two-channel model that explicitly takes into account 
the role played by the shallow bound state in the open 
channel. We show that the presence of this bound state 
implies the possibility of having a sharp minimum in 
the three-body recombination even at negative scatter- 
ing lengths, and it introduces noticeable modifications 
to the energy dependence of the low-lying trimer levels. 
As an example, we find that the trimer ground state in 
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FIG. 2; Diagrams describing: a) atom-dimer scattering, and 
b) recombination of three atoms into an atom and a dimer. 



-'^^•^Cs never hits the three-body continuum threshold at 
large values of the magnetic field, as shown in Fig. [T] 
We extract the temperature-dependent three-body re- 
combination rate for -'^'^^Cs, finding very good agreement 
with the data presented in Ref. and for ^^K, a novel 
bosonic system which appears to be a promising candi- 
date for observing Efimov effects 0] . For both species we 
also present results for the atom-dimer scattering length, 
which may be of direct relevance to future experiments. 



II. ATOM-DIMER SCATTERING 

In order to investigate three-body processes we intro- 
duce the transition matrix Ts, which describes the col- 
lision between an atom and a dimer. T3 satisfies the 
Lippmann-Schwinger equation sketched in Fig. [5^. In 
the presence of purely s-wave interactions and provided 
that the two-body T-matrix T2{E) is analytic in the 
lower complex plane, the resulting integral equation reads 

MM 



T3{k,k';E) 



1 



-Fik') 



2kk 

^^Fik")T2iE^3k"'/4m)T:,ik",k';E), 



(1) 



where E is the total energy in the center-of-mass ref- 
erence frame, k and k' are the moduli of the mo- 
menta of the incoming and outgoing atoms and F(x) — 



In [{E -k'^-x'^ + kx)/{E ~ fc2 



kx)]. The use of 



a simple T2 proportional to {1+ika)^^ in Eq. ^ presents 
a two-fold problem. In first instance, this form of the 
T-matrix does not carry any information on the pres- 
ence of a shallow level in the open scattering channel, 
which plays a fundamental role in many experimental re- 
alizations of Feshbach physics. In addition, Eq. ^ with 
T2 oc {I + ika)~^ has a well-known ultraviolet cut-off 
dependence [12j . In the following, we present a model 
which overcomes both problems. 

The atom-atom scattering length near resonance can 
be written as a{B) — abg(-B) [1 — AB/{B — Bo)], where 
Bo and AB are the location and the width of the 
resonance, and abg(S) is the scattering length in the 
open channel which may have a weak magnetic-field 



dependence [13[. Physically, the processes involved in 
the dimer propagator T2{E) can be divided into col- 
lisions occurring in the open channel only, and colli- 
sions that proceed through the formation of a molecule 
in the closed channel [3, [IB]- The low-energy col- 
lisions between the atoms in the open channel con- 
tribute the term Tbg(£;) = Tbg/ [1 - TbgS(£;)], where 
Tte = Tbg(O) = ^■nh^a^g/m and 'E.{E) = m^/'^y/^/^Tih^ 
[l^ . The propagator for a molecule in the closed chan- 
nel can be written as 1/[E — 5 — hYj{E)\, where 5 = 
Afi{B — Bq) is the energy detuning from resonance, and 
A/i is the magnetic-moment difference between the open 
and closed channels. The self-energy of the molecule is 
given by [H hT,{E) = g^E{E)/ [1 - n^EiE)], where 
g = y/Th^ApAB is the atom-molecule coupling, and 
g / [1 — ThgE{E)] corresponds to an atom-molecule ver- 
tex dressed by repeated scattering in the open channel. 
Finally, the T-matrix has to be multiplied by a factor 
of 2 to take into account the possibility of interchanging 
the identical incoming particles. By adding all the terms, 
one obtains [la, 



T2{E) 



1 



1 - TbgS(£;) J E-S- m:{E) 



AfiAB\ 



(2) 



E-6 



j - rbgS(i?) 



The pole structure of T2{E) quantitatively reproduces 
the low-energy Feshbach physics across the whole reso- 
nance. In particular, if Obg > there are two branches, 
which describe the avoided crossing between the shallow 
bound states in the open and closed channel (see Fig.[T]). 
If abg < a single branch is present, in accord with the 
fact that the open channel no longer supports a shallow 
bound state. We have explicitly checked that the eigen- 
states of T2 form an orthonormal set [l^ . 

In order to ensure the convergence of Eq. ([T]), at large 
energies T2 needs to fall-off at least as fast as 1/E [19j . 
In the absence of background interactions, i.e., if abg = 
0, Eq. (21) would behave as g'^/E and therefore would 
yield a convergent Eq. ([1]) , but we see that by adding the 
background term we obtain T2 {E) cx 1/ \J ~E. To restore 
the correct high-energy behaviour g"^ /E, we add in the 
denominator of T2 a term linear in energy that we write 
in the form a\ygR*k'^, where R* = fi^ /{mahgABAfi) > 
[2l| . In its final form, our two-body T-matrix reads 



T2{E) 



2TbB 



E-S j 



(3) 



TbgS(£;) +abgi?*/c2 



We assume here that only the two molecular states in- 
volved in the avoided crossing are important for the Efi- 
mov physics of interest, and that the effect of the much 
deeper lying molecular states can be neglected. Our T- 
matrix is expected to be accurate as long as kRo <ti 1, 
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FIG. 3: (Color online) Three-body recombination length p3 
for Cs: the curves represent our thermally-averaged results, 
while the symbols are the experimental points of Ref. ^3,] at 
lOnK (+), 200nK (x) and 250nK (□). 



where Rq = ^y—h?/mEo ~ i? is the length scale asso- 
ciated with the energy of the next lower lying state not 
included in the model. Since R* introduces in Eq. ([T]) an 
effective cut-off k* of order 1/i?*, this amounts to require 
R* :s> Ro- Our model is therefore ideally suited to de- 
scribe three-body physics around narrow resonances with 
a large background scattering length, and has no free pa- 
rameters. In case the condition R* ^ Rq should not 
be satisfied in the system under consideration, one may 
consider the quantity R* as a free parameter which 
plays then the role of an effective momentum cut-off. 
This single free parameter can then be adjusted to fit at 
best the available experimental data. This is to be com- 
pared to the 4 free parameters needed in the context of 
the universal one-channel model employed in Refs. [1, Q 
to interpret the data at positive and negative values of a. 



III. THREE-BODY RECOMBINATION RATE 

We calculate here the rate at which three atoms re- 
combine into an atom and a dimer with binding energy 
i?b- This process is sketched in Fig. We work in 
the center-of-mass frame of the three atoms and intro- 
duce the Jacobi coordinates p and q, in terms of which 
the total incoming kinetic energy is diagonal and given by 
E = /m+3q^ /Am. The T-matrix describing the recom- 
bination depends only on the magnitude of p and q, and 
is given in terms of the dimer and trimer T-matrices by 
r,ecb,g;-Bb) = 3VZT2{pym)T3ikf,q;E) 0. The factor 
of 3 in this expression comes from the sum of equivalent 
diagrams obtained by a cyclic permutation of the three 
incoming bosons, and kf — \jAm{E — i5b)/3. The result 
contains the square root of the residue Z of the dimer 
propagator evaluated at its binding energy, which takes 
into account the presence of one external dimer line. 



The three-body recombination rate is given by Fermi's 
Golden Rule, which here yields , 2kf iTj-ecP /37r. The 
sum over kf takes into account the fact that, when 
flbg > and B < Bq, one has to sum over the two possi- 
ble outcoming molecular states with binding energies E^ 
and E^ . The variation in time of the atomic number den- 
sity n is obtained by multiplying the rate by the number 
density of triples of particles in the gas n'^/3!. Assuming 
that after recombination both the atom and the dimer 
escape from the trap, we obtain dn/dt = —3Ln'^, where 
the rate constant L for a thermal gas is given by 

L = J2(^Z\T,{pym)T,{kf,q;E)A , (4) 

/th 

and (. . . )th denotes a thermal average over the initial 
Jacobi momenta. 

The recombination length p3 = {6mL/^/3hy^^ dis- 
plays narrow peaks at the magnetic field values at which 
a trimer state hits the continuum of unbound atoms. Our 
results for are compared with the experimental data 
of Ref. g in Fig. [SI Since for the broad cesium resonance 
R* — 0.15ao is not larger than the range R, we solve our 
model by keeping R* as a free parameter. We find that 
the best fit to the experimental data is obtained by choos- 
ing = 0.14ao. The width and the height of the peak, 
as well as both wings at large values of |a| are very well 
reproduced by our calculations at lOnK and 200nK. The 
recombination maximum at a = — 850ao lowers and shifts 
towards larger values of magnetic field as the tempera- 
ture is raised, in agreement with the fact that the Efimov 
state evolves into a triatomic resonance in the continuum. 
We also find a minimum of ps on the positive a side at 
a = 240ao, in agreement with the experiment. In addi- 
tion, we observe a minimum of ps at a negative value of 
the scattering length, a = — I50ao. This minimum has 
the same origin as the ones present for a > 0: it is due 
to destructive interference between two possible paths 
in the recombination process, the so-called Stiickelberg 
oscillations Indeed, caesium possesses a shallow 6s 
two-body bound state even for a < 0, and this state is 
precisely the E^ included in our theory (the lower blue 
dashed line in Fig. [T]). This novel non-universal predic- 
tion could be verified by a detailed measurement of the 
region — 200ao < a < 0. We have checked that this min- 
imum exists as long as Obg S> R, i.e., as long as E^ falls 
within the range relevant to Efimov physics. We also 
note that our results for the height of the large recombi- 
nation peak and the position of the minimum at a > 
are in much better agreement with experiment than the 
results of the parameter-free theory put forward by Lee 
et al. % 

One might wonder why i?^^ is so close to R* . The curi- 
ous proximity of the two values may be explained by the 
fact that in ^"^^Cs the next bound level lies deep down at 
Eo = -I25MHz [3, which implies i?o = 15ao < R. 
To extend the validity of our model when R* is not 
larger than Rq, we could in principle have included the 
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FIG. 4: (Color online) Recombination length ps/a at T — 
for at negative detunings. The UT fit (with decay 
parameter 77 = @|) stops at B — Bo = — 52G, where a 
vanishes and changes sign. Inset: ps/a for positive detunings 
(UT fit with a'j, given by our model and rj = 0.06). 
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FIG. 5: (Color online) Real part of the atom-dimer scattering 
length: + (-) refers to collisions with the dimer of binding 
energy (E^)- The shaded region at B < Bo contains a 
large number of resonances of ajj-,, each in correspondence 
with the appearance of a new Efimov state. The scales on the 
axes are obtained by plotting (B — Bo)^''^ and a]^-^. 



effect of this and other deeper-lying bound states by 
adding more poles to the two-body T-matrix T2{E). In 
Eq. ^ we have to integrate over these poles, which 
physically has two effects. The real, or principle value, 
part of the pole describes the virtual effects of the asso- 
ciated molecular state, whereas the imaginary part de- 
scribes the three-body recombination processes forming 
this molecule. The integral over the real part of the pole 
contains a large cancellation because of the antisymmet- 
ric nature of the integrand. Moreover, recombination 
processes forming deeply-bound molecular states are al- 
ways suppressed because of the strongly reduced overlap 
between the initial and final states. Together, these ob- 
servations indicate that sufficiently bound levels play a 
negligible role in the considered recombination process, 
and therefore they may partially explain the fact that 

^nt-^ ■ 

We now consider the three-body recombination for ^^K 
atoms around the resonance located at Bq — 403G. 
This resonance has a negative background scattering 
leng th flbg = — 29ao and is characterized by Sfi — 1.5/is 
[22| . Since for this resonance R* = 29ao is larger than 
Rq = 13ao, our model allows us to predict Efimov fea- 
tures without having to rely on existing experimental 
data. As shown in Fig. [H non-universal effects play an 
important role far from the resonance, in a region that 
could easily be accessible experimentally. Since no shal- 
low bound state exists for potassium at positive detun- 
ing 6, Eq. (HI may be used to evaluate ps only at (5 < 0. 
Nonetheless, our approach allows us to extract the un- 
known phase of the universal result at 5 > by evaluat- 
ing the scattering length a'^ at which one of the highly 
excited trimer levels hits the three-atom continuum. 



IV. ATOM-DIMER SCATTERING LENGTH 

We conclude the paper by discussing atom-dimer col- 
lisions. The atom-dimer scattering length cad is related 
to T3 by 



37rfi.^aAD 



(5) 



Here 2m/3 is the atom-dimer reduced mass and the right 
hand-side is multiplied by a factor of ^/Z for each of 
the two external dimer lines. The atom-dimer scattering 
length is displayed in Fig. [SJ and presents a resonance 
each time an Efimov state crosses the atom-dimer thresh- 
old. In the case of caesium, Obg > and two dimers are 
present at B < Bq. As a result, the scattering length for 
the collision between an atom and the shallow dimer with 
energy E — E^ is complex, since this dimer may decay 
into the deeper one. These resonances could in principle 
be observed by populating the molecular state (via RF- 
driven transitions, magnetic field sweeps, or three-body 
recombination) in '^^K, where the Feshbach resonance is 
located at > 0. 

We hope that our work will stimulate new experiments 
which are at present highly desirable to arrive at a com- 
plete microscopic understanding of the fascinating Efi- 
mov physics in ultracold Bose gases near a Feshbach res- 
onance. 
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